Critical Casimir force in the superfiuid phase: Effect of fluctuations 
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We have considered the critical Casimir force on a 4 He film below and above the bulk A point. 
We have explored the role of fluctuations around the mean field theory in a perturbative manner, 
and have substantially improved the mean field result of Zandi et al [Phys. Rev. E 76, 030601(R) 
(2007)]. The Casimir scaling function obtained by us approaches a universal constant (— gjp) for 
T<2.13K. 
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1. Introduction 

Recently, Garcia & Chan and Ganshin et al measured 
the Casimir force induced thinning of the liquid 4 He film 
near the bulk A point (T A = 2.1768 K) [J Q. They 
obtained a universal scaling function ($) for the critical 
Casimir force below and above the A point, and observed 
a dip minimum and a non vanishing constant tail in the i) 
below the A point. This experiment challenges our under- 
standing of the finite size effects on the films near their 
bulk critical points. On this issue, the Casimir effects on 
different critical films have been the subject of a number 
of experimental and theoretical (8T-fl6j| works within 
the last few years. 

Although the scaling function was appreciably ob- 
tained by the Monte Carlo simulations of Hucht [10j and 
Vasilyev et al [ll| yet this problem is still unsolved an- 
alytically. That the confinement of the critical fluctua- 
tions may give rise to a (classical) Casimir force was first 
proposed by Nightingale and Indekeu [l?}- Thereafter, a 
renormalization group calculation for the $ above the T\ 
was presented by Krech and Dietrich [l8|. For T < T\, 
a mean field theory with the Ginzburg-Landau (G-L) 
model was recently presented by Zandi et al [3]. They 
obtained an analytic expression for the $ in terms of the 
maximum of the superfiuid order parameter. By propos- 
ing that their mean field calculation could be improved 
by the confinement of the critical fluctuations (at the 
Gaussian level) , they nicely improved their result only at 
the A point. 

We analytically improve the mean field result of Zandi 
et al as proposed by them for T <T\ [9] . The improve- 
ment for T > T\, was already done by Krech and Diet- 
rich even beyond the Gaussian level [18j . However, we 
present a physically motivated regularization technique 



for obtaining the critical Casimir force above the A point. 
Thus we build a unified picture for the theory of critical 
Casimir force acting on a 4 He film below and above the A 
point. Our theory interestingly predicts the non vanish- 
ing constant tail of ■& as —0.0478, which agrees well with 
the numerical result of Hucht [10j but differs by a factor 
of five from the experimental value (—0.24) Q. Nonethe- 
less, it is a considerable improvement over the mean field 
calculation which predicts it to be zero @ . 



We start from the G-L model. For T > T\, we ob- 
tain the free energy in terms of the discrete Fourier 
modes. The Casimir force is then obtained in the Fisher- 
de Gennes' form [l9| by applying the Poisson summa- 
tion formula [2oT ]. Use of this summation formula dis- 
tinguishes our approach from that of Krech and Dietrich 
[18] . For T <T\, we transform the critical fields by intro- 
ducing the superfiuid order parameter, and express the 
G-L free energy in a decoupled form of the mean field 
and fluctuating parts. The fluctuating part is treated 
like that we do for T > T\, and the mean field part is 
treated in the manner of Zandi et al Q . It is necessary to 
know the maximum of the order parameter for plotting 
the mean field part of the Casimir force. Although the 
graphical solutions of the maximum of the order param- 
eter are exact yet the solutions do not appear in a closed 
form. We predict a closed form of the maximum of the 
order parameter from asymptotic analyses, and obtain 
an approximate mean field Casimir force which matches 
very well with the exact mean field result 9] . Finally we 
improve the mean field result by adding the contribution 
of the fluctuating part. 
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2. Free energy of the critical fluctuations for T > T\ 

According to the experimental setup 4 He vapor comes 
in contact of a plate, and upon liquefaction it forms a 
film of thickness 238-340A 0,0]. We consider the plate 
to be along the x — y plane of the co-ordinate system, the 
area of the him to be A, and the thickness of the film to 
be L along the z direction. Near the A point 4 He behaves 
critically, and its local free energy can be written in the 
G-L framework as 
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where <f>(r) = <fii(r) + i(fc(r) is a complex scalar critical 
field at the position vector r = xi+yj+zk, a is the inverse 
square of the correlation length (£ = £ot~ u ), t = T/T\ — l 
is the reduced temperature, v is the correlation length 
exponent, and & is a positive coupling constant pH liij]. 
The quartic term in Eqn.(l) is neglected in the Gaussian 
approximation. 

Let us first calculate the Casimir force for T > T\, 
In conformity with the Dirichlet boundary conditions, 
the Fourier expansion of the critical fields are given by 

*i,2to = yf Er=i/^,2,»(k)sii 1 (^)e ik -(-+«) I ^. 
In the basis of the Fourier modes, we obtain the par- 
tition function (Z = j D[4> 1 ]D[(j)2]e^ Fl / kBT ) within the 
Gaussian approximation, and get the standard form of 
the free energy (— ksTlnZ) of the critical fluctuations of 
the film as uM 



F = 2xM^f; Tlnffc 
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The factor 2 in the above equation comes from the fact 
that (j) has two components. 



where T\ — 35- ■ Using the Poisson summation formula 
in Eqn.(4), we recast Eqn.(3) as 
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As L oo, only the first term of the parentheses 
of Eqn.(5) survives. This is the bulk force acting on 
the film. By the standard analytic continuation tech- 
nique we get the expression of this bulk force as /oo = 
2x ) 3/2 r(-3/2). Subtracting this bulk part from 

fr, we get the Casimir force in the Fisher-de Gennes form 
19] f c [L,t] = ^0±ti{t) where §{t) is the Casimir scal- 
ing function which can be expressed in terms of a scaled 
temperature r = L x l v t as 
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We have v = 1/2 in the Gaussian as well as in the mean 
field approximations [22[. However, if we want to in- 
clude the effect of the |(/>| 4 term within the above pre- 
scription, we must put v = 0.67016 (~ 2/3) in Eqn.(6) 



[E l2l, |22|. From Eqn.(6) we get the value of the $(r) at 

the A point as — = —0.0956 which matches well with 
• • • l~Ti 

the experimental data obtained by Garcia and Chan [1|. 

The same number at the A point was also obtained in 

Refs.fil. [l8l| with different regularization techniques. We 

need to know the value of £o for plotting the $(t) against 

t. The experimental value of £o for T > T\, varies from 

1.2 to 1.43A [13-ii]. With no a priory reason we take 

Co = 1-3A for T>T X [24|. 



4. Free energy of the critical fluctuations for T <T\ 



3. Critical Casimir force for T > T x 

From Eqn.(2) we get the force acting on the film as 

_ DF _ ATA 

~~ ~dL ~ L 3 (3) 

where S = J2n°=i Iq° ~ — ""a^ This expression can 
be recast as 

ikJo Jo ^ 



In addressing the situation below the A point we note 
that the Casimir scaling function looks qualitatively sim- 
ilar to the ultrasonic attenuation and finite size specific 
heat (i.e. both have a peak below T A ) [H, H3]. We 
anticipate that the Casimir effect for T < T\ can be 
thought of as coming from mean field and fluctuating 
parts. Splitting the ultrasonic attenuation into a sum of 
mean field and fluctuating parts was the original contri- 
bution of Landau and Khalatnikov, and it gives a good 
account of the ultrasonic attenuation below the A point 
[28| . Here we show how a similar approach for the critical 
Casimir effect can be adopted below the A point. 

We return to Eqn.(l) and note that for T < T\, a 
becomes negative, and accordingly we write a = —\a\. 
This leads to a broken-symmetry, and we handle it by 
transforming the fields <f>2 and 4>i to tp2 = <\>2 and 
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ipi = (f>i—m{z) where m{z) is the superfluid order param- 
eter. The expectation value < <f>\ > is now z-dependent 
because we are considering a finite size system in the z- 
direction and consequently, we expect an inhomogeneity 
in the superfluid density (~ m 2 (z)). The fields ipi,ip2 are 
such, that < tfji >— (i = 1, 2), and the local free energy 
in Eqn.(l) in terms of V'IjV^ becomes 
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The free energy can be minimized from the condition that 
< >= 0, and can be recast from Eqn.(7) as [— ^£ — 
\a\m + bm 3 } +bm[< ip% > +3 < ip 2 >] = which can only 
be solved analytically if we disregard the second square 
bracketed term by considering the necessary condition 
that the mean field part dominates over the fluctuating 
part (i.e. m 2 ^< ipf >). With this consideration we can 
write an approximate equation for the profile of m(z) as 
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(8) 



It is to be noted that Eqn.(8) does not minimize the lo- 
cal free energy in Eqn.(7). Hence, the quadratic terms 
in ipi and ip2 may not be positive. However, Eqn.(8) 
would minimize the local free energy if we replace m(z) 
in the quadratic and higher order terms by its bulk value 
(i/|a|/6). With all the above considerations (and with 
Eqn.(8)) the fluctuating and mean field parts of the lo- 
cal free energy become decoupled, and consequently, we 
recast Eqn.(7) as 
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is the 
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where F m j 
mean field part, F a 

^(V-02) 2 ] is the (Gaussian) fluctuating part, and F, ? 
Jd 3 (r)[ N /ja|6?/'i(V'i + ip%) + |(^ 2 + ^f) 2 ] is the interac- 
tion part. We now check that all the quadratic terms in 
Fi are positive. Hence, Eqn.(8) minimizes the local free 
energy, and consequently, the quadratic terms in Eqn.(9) 
on an average dominate over the higher order terms as 
because ^ ~ m 2 ^$>< ip? >. Evaluation of the partition 
function from the extremized local free energy in Eqn.(9) 
leads to 
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where Z = J D[ipi]D[ip2\& kBT is the partition function 
for the (Gaussian) fluctuating part, and the expectation 
value < ... > is taken with respect to F . The free 
energy obtained from Eqn.(10) is given by 
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where F c f = — ksT\\vZ is the (Gaussian) fluctuating 
part of the free energy for T <T\. 



5. Critical Casimir force for T <T\ 

5. A. Fluctuating contribution 

The fluctuating part of the free energy in Eqn.(ll) 
can be recast in a special form of Eqn.(2) as F c f — 

EZi Io°° [H? + 2\a\ + 4^] + \n[k 2 + 4^]] 4f 
which gives the Casimir scaling function 
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by following the steps from Eqn.(2) to Eqn.(6). 

Since ^ s» m 2 ^>< ipf > we can easily check from 
Eqn.(ll) that F mf > F cf >< F int > Q . Thus we can 
ignore the < Fi nt > terms in Eqn.(ll), and can expect 
that the Casimir force obtained from the fluctuating part 
would be much smaller than that obtained from the mean 
field part. 



5.B. Mean field contribution 

Let us now evaluate the Casimir scaling function from 
the mean field part (F m f) in Eqn.(ll). From the con- 
sideration that the order parameter m(z) is smooth and 
obey the Dirichlet boundary conditions m(0) = m(L) = 
0, m(z) must be symmetric about z = h and there would 
be a single maximum of m(z) at z — U for the lowest pos- 
sible value of the mean field free energy. An analytical 
expression of the Casimir force (— [ — — ~gj^ Il->oo] ) 
from the mean field part F m f (and from Eqn.(8)) was 



nicely obtained in Ref. [9j in terms of rj 
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FIG. 1: The dots are the graphical solutions of 77 in Eqn.(14). 
The dotted and dashed lines are the asymptotic solutions near 
L\/Ja\ — > 7r and L^/\a\ 00 respectively. The continuous 
line follow Eqn.(15). 



In Eqn.(13) rj as well as the maximum of the order pa- 
rameter is restricted by 0, [2^] 



length as 



fmf — 
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for 7T > Lji > 



for Lji > ir. 



(16) 

From Eqn.(16) a dip minimum with discontinuous slope 
is expected to occur at Ly/\a~\ = ir. This point is fitted to 
the experimental dip at r = —9.7 2]. The modifications 
to Eqn.(16) would come from the higher order fluctuat- 
ing terms, and the primary correction would be to keep 
the form of the /,„/ unaltered with the mean field £ be 
replaced by £ = £ol^l~ 1 ' where v to the lowest order in b 
is i + |(n + 2) [22]. Using the fixed point value of b we 
can get the usual v at one loop order. We can safely as- 
sume that the effect of the different loops will be to make 
i = i Q t- v with v acquiring the value 0.67016 (« 2/3) 
correct to all orders [2, 111]. From Eqn.(16) we get the 
modified mean field Casimir scaling function in terms of 
t = L x l v t as 



L 



( 14 ) t? m/ (r) = 



64 X% 
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where K(x) 

plete elliptic integral of the first kind, 
the allowed range < rj < \ for the corresponding do- 



- ...] is the com- 
Eqn.(14) gives 



main n < L^/H < 00, and it can be exactly solved 
by the graphical method (29^ . Although the graphical 
method does not provide 77 in a closed form of L-y/jaj yet 
we can do so by the asymptotic analyses near 77 — > 
and r\ — > \. For L-\/[a| — > 00, the asymptotic so- 
lution of r) in Eqn.(14) is r) -> itanh 2 (^p^) 0. 
On the other hand, for L-v/|a| — > it, the asymptotic 



solution (up to the third order in rj in Eqn.(14)) is 
^-^§(W(l-i(^) 2 + l-04514(^) 4 +...) where 
a = |a| — J? (29| |. Corresponding to the above asymptotic 
solutions we can take a fitting function for the domain 
7T < Ly^T\ < 00 as [H, [H 



, r /TTV 1 u2 ^(£ 2 | a |-^ 2 )/2 

r/(iVl a l) = 2 tanh [ 2 



(15) 



We see in FIG. 1 that all the asymptotic and graphical 
solutions match very well with Eqn.(15). Hence, we con- 
sider Eqn.(15) as an approximate solution for the rest of 
this paper. 

With the consideration of Eqn.(15), and that 77 = [9( 
for < Ly^jaj < 7r, one can recast Eqn.(13) in terms of 
the reduced temperature and the mean field correlation 



H 2 



for - 9.7 < t < 



s.vli I \/ — ^ I for r < -9.7. 

(17) 

If we plot Eqn.(17) we must get almost the same result as 
obtained by Zandi et al [ij. However, we need to improve 
the mean field result (in Eqn.(17)) by the confinement of 
the critical fluctuations as proposed by them. 



5. C. Improvement of the mean field result 

The net Casimir scaling function for T < T\ is ob- 
tained from Eqns.(12) and (17), and is given by 



0(T) = m /(r) + i?o/O-)- 



(18) 



We plot the right hand sides of Eqns.(6) and (18) in FIG. 
2. For T > T\, our theory matches very well with the 
experimental data of Garcia & Chan [lj. From FIG. 2 
we also see that the $ approaches a constant 



03) 
8tt 



-0.0478 for r < -75.6 A ' and (with L = 238A) for 
T < 2.13 K as well. 

Although our theory for T < T\, does not match very 
well with the experimental data yet it predicts the basic 
nature of the critical Casimir force which is characterized 
by a non vanishing constant tail jH • It is of course clear 
from FIG. 2 that the inclusion of the effect of the critical 
fluctuations substantially improves the exact mean field 
result of Zandi et al [9j . 
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FIG. 2: The continuous line for r > 0, follows Eqn.(6) with 
fo = 1-3A and that for r < 0, follows Eqn.(18) with £ = 
1.4593A, AKci } = = 0.0133547 and v = 0.67016. The dashed 

line represents exact mean field calculation of Zandi et al [y| 
with Co = 1.4593A, , Kci ) = = 0.0138166 and v = 0.67016. 

A few experimental points are taken from Ref. [1] for r > 0, 
and from Ref.Q] for r < 0. 



The experimental dip of the # has been adjusted with 
the value of b which nobody has determined (for the film) 
so far from the theoretical point of view. How to deter- 
mine the parameter b for the film and to calculate the 
Casimir force by considering the coupling between the 
mean field and fluctuating parts remain to this day as 
open problems. 
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6. Conclusions 



Complemen ting the numerical 1(], II, 14 1 and analyt- 
ical @, IH E, EI El, El works we have given a unified 
theory for the the critical Casimir force below and above 
the lambda point in a single framework. In particular, 
we have explored the effect of the critical fluctuations in 
the Gaussian level over the mean field contribution 

The tail of the scaling function approaches a non- 
vanishing constant — = —0.0478 owing to the con- 
sideration of the two component (ipi,ip2) critical fluctu- 
ations. Although this constant is closer to the numerical 
simulation result obtained by Hucht [10( yet it is nearly 
one fifth of that obtained by the experimentalists [2j . On 
the other hand, this constant is zero in the mean field 
level [9(. Hence, our calculation of the Casimir scaling 
function goes beyond that of Zandi et al [9[, and com- 
pares favorably with the numerical simulation of Hucht 
and the experimental data of Garcia & Chan [l| and 
Ganshin et al [2j. 

It should be mentioned that the theory of the Casimir 
force for T < T\, was also improved by Zandi, Rudnick 
and Kardar with the considerations of confinement of the 
Goldstone modes & surface fluctuations [3(| , and that the 
thinning of the liquid 4 He film was first (but admitt edly 
not very precisely) observed by Dionne and Hallock [31| . 

While the Casimir force for the quantum fluctuations 
of the electromagnetic field is observed within 10 -12 N 
[32| the critical Casimir force considered by us is ob- 
served within 10 -3 N [l], 0]. The confinement of the 
classical (critical) fluctuations of course is much stronger 
than that of the quantum (vacuum) fluctuations. 
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